Effect of large-scale coherent structures on turbulent convection 
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We study an effect of large-scale coherent structures on global properties of turbulent convection 
in laboratory experiments in air flow in a rectangular chamber with aspect ratios A~ 2 and A « 4 
(with the Rayleigh numbers varying in the range from 5 x 10^ to 10*). The large-scale coherent 
structures comprise the one-cell and two-cell flow patterns. We found that a main contribution to 
the turbulence kinetic energy production in turbulent convection with large-scale coherent structures 
is due to the non-uniform large-scale motions. Turbulence in large Rayleigh number convection with 
coherent structures is produced by shear, rather than by buoyancy. We determined the scalings of 
global parameters (e.g., the production and dissipation of turbulent kinetic energy, the turbulent 
velocity and integral turbulent scale, the large-scale shear, etc.) of turbulent convection versus the 
temperature difference between the bottom and the top walls of the chamber. These scalings are in 
an agreement with our theoretical predictions. We demonstrated that the degree of inhomogeneity 
of the turbulent convection with large-scale coherent structures is small. 

PACS numbers; 47.27.te, 47.27.-i 



I. INTRODUCTION 

Coherent structures in a turbulent convection at large 
Rayleigh numbers have been observed in the atmospheric 
turbulent convection [H, i, S i, i, i, 0, S i, 0, [lH E 

■ [l3 | (so-called "semi-organized structures" ), in numer- 
ous laboratory experiments in the Rayleigli-Bcnard ap- 
paratus |23|, IM y, [l, 

[23, |2i, [2a, [30, [Ml, IM [33, m, m, m, (so-called "wind" 
or "mean wind") and in direct numerical simulations 
[stI . [bsI . [sol . [40l |. Characteristic spatial and time scales 
of the coherent structures in a turbulent convection are 
larger than turbulent scales. 

In spite of numerous theoretical, experimental and nu- 
merical studies of coherent structures in turbulent con- 
vection, their origin is still a subject of numerous discus- 
sions. In particular, there are two points of view on the 
origin of coherent structures in turbulent convection [ssl ] . 
According to one point of view, the rolls which develop 
at low Rayleigh numbers near the onset of convection 
continually increase their size as Rayleigh number is in- 
creased and continue to exist in an average sense at even 
the highest Rayleigh numbers reached in the experiments 
(4l| . Another opinion consists in that coherent structures 
are genuine high Rayleigh number effect (Tsj . 

A new mean-field theory of formation of coherent 
structures in turbulent convection has been suggested 
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in [ij, [43, [43|. According to this theory a redistribu- 
tion of the turbulent heat flux due to non-uniform large- 
scale motions plays a crucial role in the formation of the 
large-scale coherent structures in turbulent convection. 
A convective-wind instability in the shear-free turbulent 
convection causes the large-scale motions in the form of 
cells. In the sheared convection, the large-scale instabil- 
ity results in an excitation of convective-shear waves, and 
the dominant coherent structures in the sheared convec- 
tion are rolls [ll[4l,[4l|. 

It was suggested on the base of laboratory experiments 
and numerical simulations performed in [23 that the co- 
herent structures in the Rayleigh-Benard turbulent con- 
vection are not driven by the turbulent Reynolds stresses, 
associated with the tilting plumes at the upper and the 
lower horizontal walls. The numerical results in [23| show 
that once the mean flow is established, the temperature 
of the fluid is larger at one side wall and smaller at the 
other side, and the mean flow is driven by the mean buoy- 
ant force at the side walls. This is also in agreement with 
the experimental studies in (25l . [soj . 

There is an opinion that buoyancy plays a crucial role 
in production of turbulent convection. For example, in 
the mixing length theory of astrophysical turbulent con- 
vection the followin g es timate for the turbulent kinetic 
energy is often used [il. [isl [4^: 



{xi^)/2=gT{u,s), 



(1) 



where {uzs) is the vertical turbulent flux of entropy, u 
and s are fluctuations of fluid velocity and entropy, g is 
the acceleration of gravity, r is the characteristic correla- 
tion time of turbulent velocity field and angular brackets 
denote ensemble averaging. Equation ([T]) implies that 
the vertical turbulent flux of entropy plays a role of a 
stirring force for the turbulent convection. 

However, this estimate can be valid only in the ab- 
sence of large-scale coherent structures. Measurements 
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of the mean temperature distributions and heat fluxes 
in laboratory experiments in the Rayleigh-Benard tur- 
bulent convection (see, e.g., [H, [sil, l47( ) show that a 
thermal structure inside the large-scale circulation is in- 
homogeneous and anisotropic. The hot thermal plumes 
accumulate at one side of the large-scale circulation, and 
cold plumes concentrate at the opposite side of the large- 
scale circulation. These measurements demonstrate that 
the vertical turbulent heat flux inside the large-scale cir- 
culation is very small in spite of large temperature dif- 
ference AT between the bottom and the top walls of the 
chamber. On the other hand, the horizontal turbulent 
heat flux inside the large-scale circulation is larger than 
the vertical turbulent heat flux in spite of the absence of 
the imposed temperature difference (AT)j, between the 
side walls of the chamber (see, e.g., |28i]). The mean ve- 
locity field inside the large-scale circulation is strongly 
non-uniform and turbulence can be also produced by 
mean velocity shear. Therefore, one may ask the question 
about the origin of turbulence inside the coherent struc- 
tures, i.e., is it shear-produced turbulence or buoyancy- 
produced turbulence? 

The goal of this study is to investigate experimentally 
an effect of large-scale coherent structures on global prop- 
erties of turbulent convection. In particular, we address 
the following issues: (i) the origin of production in tur- 
bulent convection with large-scale coherent structures (a 
shear-produced or buoyancy-produced turbulence); (ii) 
the scalings of global parameters (production and dissi- 
pation of turbulent kinetic energy, turbulent velocity and 
integral turbulent scale, the large-scale shear, etc.) ver- 
sus the temperature difference AT between the bottom 
and the top walls of the chamber; (iii) the degree of in- 
homogeneity of the turbulent convection with large-scale 
coherent structures. 

The paper is organized as follows. Section II describes 
the experimental set-up for a laboratory study of the co- 
herent structures. The experimental results and their de- 
tailed analysis are presented in Section III. Finally, con- 
clusions are drawn in Section IV. 



II. EXPERIMENTAL SET-UP 

The experiments were conducted in two rectangular 
chambers with dimensions 26 x 58 x 26 cm'^ and 26 x 
58 X 13 cm'^. Hereafter, we use the following system of 
coordinates: z is the vertical axis, the y-axis is directed 
along the longest wall (see Fig. [T]). The side walls of 
the chambers are made of transparent Perspex with the 
thickness of 10 mm. A number of experiments have been 
conducted with different additional thermal insulation of 
the side walls of the chamber in order to study whether 
a heat flux through the side walls affects the turbulent 
convective pattern. First, the side walls of the chamber 
were insulated with Styrofoam plates with low thermal 
conductivity (k ~ 0.033 W/mK) and with the thickness 
of 30 mm. Two of the side plates were removed for a 




FIG. 1: Experimental set-up: (1) - walls of the chamber; (2) 
- cooled top with a heat exchanger; (3) - laser light sheets; 
(4) - heated bottom with a heat exchanger; (5) - two digital 
CCD cameras. 



short time when the images of the flow were recorded. 
In the next series of experiments we installed additional 
Perspex plates with a thickness of 6 mm and with an 
air gap of 2 mm between these plates and the outside 
walls of the chamber. Finally, we performed experiments 
where these two types of thermal insulation were used 
simultaneously. All these experiments have shown that 
the heat flux through the side walls does not affect the 
turbulent convective pattern. 

A vertical mean temperature gradient in the turbulent 
air flow was formed by attaching two aluminium heat 
exchangers to the bottom and top walls of the test section 
(a heated bottom and a cooled top wall of the chamber) . 
The thickness of the aluminium plates is 2.5 cm. The 
top plate is a bottom wall of the tank with cooling water. 
Temperature of water circulating through the tank and 
the chiller is kept constant within 0.1 K. Cold water is 
pumped into the chamber through two inlets and flows 
out through two outlets located at the side wall of the 
chamber. The bottom plate is attached to the electrical 
heater that provides constant and uniform heating. The 
voltage of a stable power supply applied to the heater 
varies up to 200 V. The power of the heater varies up 
to 300 W. The temperatures of the conducting plates 
were measured with four thermocouples attached at the 
surface of each plate. 

The temperature difference between the top and bot- 
tom plates, AT, varies in the range from 20 K to 80 K 
depending on the power of the heater (i.e., the global 
Rayleigh number based on molecular transport coeffi- 
cients, was changed within the range of Ra = (0.4 — 1) x 
10® for the aspect ratio of the chamber A ~ Hy/Hz — 
2.23 and Ra = (0.5 - 1.2) x 10^ for A = 4.46). Here 
Hy and Hz are the sizes of the chamber along y-axis 
and z-axis, respectively. The temperature in the probed 
region was measured with a specially designed temper- 
ature probe with twelve sensitive thermocouples having 
the diameter 0.125 mm. This allows us to obtain the de- 
tailed mean temperature distributions inside the coherent 
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structures. 

The velocity field was measured using a digital Particle 
Image Velocimetry (PIV) , see [i^ . In the experiments we 
used LaVision Flow Master III system. A double-pulsed 
light sheet was provided by a Nd-YAG laser (Contin- 
uum Surelite 2 x 170 mJ). The light sheet optics includes 
spherical and cylindrical Galilei telescopes with tuneable 
divergence and adjustable focus length. We used two 
progressive-scan 12 bit digital CCD cameras (with pixel 
size 7.4 ^m x 7.4 ^m and 2048 x 2048 pixels) with a dual- 
frame-technique for cross-correlation processing of cap- 
tured images. A programmable Timing Unit (PC inter- 
face card) generated sequences of pulses to control the 
laser, camera and data acquisition rate. The software 
package LaVision DaVis 7 was applied to control all hard- 
ware components and for 32 bit image acquisition and 
visualization. This software package comprises PIV soft- 
ware for calculating the flow fields using cross-correlation 
analysis. The velocity field was measured in the two per- 
pendicular planes (in five yz planes and in ten xz planes, 
so that the distance between these planes is 5 cm, see 
Fig.[l]). 

An incense smoke with sub- micron particles as a tracer 
was used for the PIV measurements. Smoke was pro- 
duced by high temperature sublimation of solid incense 
particles. Analysis of smoke particles using a microscope 
(Nikon, Epiphot with an amplification of 560) and a PM- 
300 portable laser particulate analyzer showed that these 
particles have an approximately spherical shape and that 
their mean diameter is of the order of 0.7/xm. The proba- 
bility density function of the particle size measured with 
the PM300 particulate analyzer was independent of the 
location in the fiow for incense particle size of 0.5 — 1 /im. 

Series of 260 pairs of images acquired with a frequency 
of 1 Hz were stored for calculating the velocity maps and 
for ensemble and spatial averaging of turbulence charac- 
teristics. The center of the probed fiow region coincides 
with the center of the chamber. We measured the velocity 
field in different fiow areas, e.g., in the chamber with the 
aspect ratio A w 2, the probed fiow area is 492 x 212.5 
mm^ with a spatial resolution of 302 /xm / pixel in yz 
plane, and the probed flow area is 237 x 211 mm^ with 
a spatial resolution of 264 /im / pixel in xz plane. Sim- 
ilarly, in the chamber with the aspect ratio A « 4, the 
probed flow area is 547 x 127 mm^ with a spatial res- 
olution of 278 /im / pixel in yz plane, and the probed 
fiow area is 256 x 125 mm^ with a spatial resolution of 
130 /im / pixel in xz plane. These regions were analyzed 
with interrogation window of 32 x 32 pixels in the cham- 
ber with the aspect ratio A k, 2 and with interrogation 
window of 24 x 24 pixels in the chamber with the aspect 
ratio A « 4. A velocity vector was determined in every 
interrogation window, allowing us to construct a velocity 
map comprising up to 3655 vectors. 

Mean and r.m.s. velocities, two-point correlation func- 
tions and integral scales of turbulence were determined 
from the measured velocity fields. The mean and r.m.s. 
velocities for each point of the velocity map (up to 3280 



points) were determined by averaging over 260 indepen- 
dent maps, and then over 3280 points. The two-point lon- 
gitudinal correlation functions of the velocity field were 
determined for the central part of the velocity map by 
averaging over 260 independent velocity maps. In the 
experiments we evaluated the variability between the fist 
and the last 20 frames. Since no variability was found, 
these tests showed that 260 image pairs contain enough 
data to obtain reliable statistical estimates. The charac- 
teristic turbulence time in the experiments Tj, = 0.8 — 1.2 
seconds, while the characteristic time period for the 
large-scale circulatory flow is by one order of magnitude 
larger than r^. These two characteristic times are much 
smaller than the time during which the velocity fields are 
measured. 

An integral scale i of turbulence was determined from 
the two-point correlation functions of the velocity field. 
These measurements were repeated for various temper- 
ature differences between the bottom and the top walls 
of the chamber. The size of the probed region did not 
affect our results. Similar experimental set-up and data 
processing procedure were used in experimental study 
of hysteresis phenomenon in turbulent convection in [3l| 
and in 49, 50, 51, 52] for investigating a new phenomenon 
of turbulent thermal diffusion |53l . l54l | . 

The maximum tracer particle displacement in the ex- 
periment was of the order of 1/4 of the interrogation 
window. The average displacement of tracer particles 
was of the order of 2.5 pixels. The average accuracy of 
the velocity measurements was of the order of 4% for 
the accuracy of the correlation peak detection in the in- 
terrogation window of the order of 0.1 pixels (see, e.g., 

ISlHElli). 

In order to measure the vertical turbulent heat ffux 
Fz = {uzO), the velocity Uz and temperature 9 ffuctu- 
ations are determined using the temperature measure- 
ments by a thermocouple simultaneously with the veloc- 
ity measurements by the PIV system. A time constant of 
the thermocouple was evaluated experimentally as 75 ms 
that is close to the Kolmogorov's time of turbulence. The 
temperature is acquired in a time interval between two 
laser pulses used for the acquisition of every velocity vec- 
tor map. During this time interval we obtain the data set 
consisting of about 100 measured temperatures. In order 
to decrease a noise level we use the mean value of this 
data set. Simultaneous temperature and velocity mea- 
surements should be conducted at the same point in order 
to determine the vertical turbulent heat flux Fz- Since 
the distance between the light sheet and the thermocou- 
ple is 3 mm, we have to use a correction factor 1/0.95 to 
a measured value of the heat flux. This correction fac- 
tor is obtained using the measured two-point correlation 
function of the velocity field. Indeed, the magnitude of 
the longitudinal correlation function of the velocity field 
at the distance 3 mm, is 0.95. 
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III. PRODUCTION AND DISSIPATION IN 
TURBULENT CONVECTION 

Our velocity measurements show that the observed 
large-scale coherent structures comprise the one-cell and 
two-cell flow patterns. For instance, in Figs. [Hand [3] we 
show the mean flow patterns obtained in the experiments 
with turbulent convection in the chambers with A « 2 
and yl « 4. The prevailing observed flow pattern in the 
experiments in the chamber with the aspect ratio A ^ 2 
is one-cell flow structure (see Fig. [5]), while the prevailing 
observed flow pattern in the experiments in the chamber 
with the aspect ratio ^ « 4 is two-cell flow structure 
(see Fig. [3]). The prevailing observed flow pattern is the 
flow pattern which is observed in the largest range of the 
Rayleigh numbers. This feature is in agreement with the 
theoretical predictions made in [43l |. In particular, the 
threshold required for the excitation of the large-scale 
instability that causes formation of large-scale coherent 
structures, is minimum when the ratio of the horizontal 
size to the vertical size of the large-scale circulation is 
approximately 2 [i^ . This value corresponds to the one- 
cell flow pattern in the experiments with A ^ 2, and to 
the two-cell flow pattern in the experiments with A ^ A. 
Note also that depending on the temperature difference 
AT between the bottom and the top walls of the cham- 
ber, we observe in the chamber with A « 4 the two-cell 
flow patterns with the downward (for AT < 26 K) or 
upward (for AT > 26 K) motions in the central region of 
the chamber between two cells (see Fig. [3]). 

Note that the large-scale circulations in a turbulent 
convection have been previously experimentally studied 
in [s^l in a Rayleigh-Benard cell with the aspect ratio 
A = 4 (with sizes 60 x 60 x 15.5 cm^^). Experiments con- 
ducted in this cell filled with water (with Prandtl number 
Pr= 5.5) and at a global Rayleigh number Ra= 5.9 x 10* 
(corresponding to a temperature difference of 6.3 K), 
revealed that the mean flow fleld measured with PIV 
comprises two rolls. This result was confirmed in our 
experiments for the similar aspect ratio of the chamber 
A = 4.45, but in air flow (with Prandtl number Pr= 0.7) 
and with different horizontal cross section of the cell: 
rectangle in our experiments and square in the experi- 
ments reported in [30|. 

The mean temperature measurements show that the 
thermal structure inside the large-scale circulation is in- 
homogeneous and anisotropic. For instance, in Fig. [4] 
we show the isolines of the mean temperature field ob- 
tained in the experiments with turbulent convection in 
the chamber with ^ w 4 in the yz plane for the two-cell 
flow patterns with the upward motions in the central re- 
gion of the chamber between two cells. The hot thermal 
plumes accumulate at one side of the large-scale circula- 
tions (in the central part of the flow), and cold plumes 
concentrate at the opposite side of the large-scale circu- 
lations (in the periphery part of the flow near side walls) . 
This fact is in agreement with the temperature measure- 
ments performed in (28j . 
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FIG. 2: Mean flow patterns obtained in the experiments with 
turbulent convection in the chamber with A ~ 2: one-cell 
flow patterns with the counter clockwise motions at the tem- 
perature difference AT = 20 K (upper panel); one-cell flow 
patterns with the clockwise motions at the temperature dif- 
ference AT = 35 K (lower panel). Coordinates y and z are 
measured in mm. 
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FIG. 3: Mean flow patterns obtained in the experiments with 
turbulent convection in the chamber with A « 4: two-cell flow 
patterns with the downward motions in the central region of 
the chamber between two cells at the temperature difference 
AT = 20 K (upper panel); two-cell flow patterns with the 
upward motions in the central region of the chamber between 
two cells at the temperature difference AT = 33 K (lower 
panel). Coordinates y and z are measured in mm. 

Now let us discuss the origin of the production and dis- 
sipation in turbulent convection with large-scale coherent 
structures. Equation for the evolution of the turbulent 
kinetic energy Ek = (u^)/2 reads: 

—^ + div^K = -{u,u,)V,U^+l3F,-e , (2) 
where D/Dt — + U-V, u are the fluctuation of the 
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FIG. 4: Mean temperature field obtained in the experiments 
with turbulent convection in the chamber with A ~ 4 in the yz 
plane for the two-cell flow patterns with the upward motions 
in the central region of the chamber between two cells. The 
mean temperature field is shown in degrees of C. Coordinates 
y and z are measured in mm. 

fluid velocity, U is the mean velocity that describes co- 
herent structure, is the term that includes the third- 
order moments, (3 = g /T* is the buoyancy parameter and 

is a reference value of the mean absolute temperature, 
g is the acceleration of gravity, Fi = {uiO) is the heat flux, 
9 are the temperature fluctuations, e = C^EkIt is the 
dissipation rate of the turbulent kinetic energy. While 
the magnitude of the constant Cg, based on several ex- 
periments that were summarized in (STj and on the data 
obtained in [ssl ] , exhibits a relatively large scatter in the 
range of Reynolds numbers Rca < 50, on the average 
the value of the constant is closed to 1. Here Re a 
is the Reynolds number based on the Taylor microscale. 
The value of the constant = 1 was obtained for the 
first time in (soj based on the experimental study 60] for 
14 < Rca < 41. 

In order to determine experimentally the turbulent ki- 
netic energy dissipation rate e, we take into account a 
small anisotropy of the turbulent velocity field. In par- 
ticular, we determine e as follows: 

e=9±l<^±^±\^ (3) 

2 V ^^z ' 

where u^^ Uy and Uz are the components of the turbulent 
velocity field (rms), £x, (y and iz are the integral turbu- 
lent length scales along x, y and z axis, and Ce = 1. In 
our experiments we determine the turbulent length scales 
and three components of turbulent velocity fluctuations 
in the range of Rca from 30 to 55, and for the temper- 
ature difference range AT between the bottom and the 
top walls of the chamber from 19 K to 80 K. We mea- 
sure components of turbulent velocity fluctuations field 
in 15 planes as indicated above, and for each velocity 
field we determined the average values of over the 
plane. For calculating integral turbulent length scales 
£i we determine the longitudinal correlation functions of 
the turbulent velocity field in each plane. These correla- 
tion functions are approximated by exponential function 
and integral scales are calculated by integration of these 
functions. 

It must be noted that the turbulent kinetic energy dis- 
sipation rate is different for various components of tur- 
bulent velocity fluctuations although the difference is not 



large. Indeed, the contribution to e of the y-component 
is 0.27 while contributions of the turbulent velocity com- 
ponents in the xz plane are 0.35 and 0.38 for x and z 
components, respectively. When turbulent kinetic energy 
dissipation rate is calculated through a one-dimensional 
surrogate Uy, as it is done in most experimental investi- 
gations, then the magnitude of the constant should 
be multiplied by a factor 1.23 which is the ratio of 1/3 
and 0.27. The magnitude of the renormalized constant 
Ce is very close to that obtained in [sH. It is conceiv- 
able to suggest that one of the reasons of the scatter in 
the measured value of the constant in the cited above 
studies is caused by anisotropy of the turbulent velocity 
fleld. 

Using the data obtained from the measured velocity 
fleld in many planes we determine the shear-induced pro- 
duction term P = — (wj Uj) V jUi of the turbulent kinetic 
energy. The main contributions to the production term 
are due the following terms: 

P = V^^^ {{VzVyf + (V,C/,)2] + V^^ [V^Uzf 

+ V^^{VyUzf + .... (4) 

where v^^ = Uxix and similarly for v^^^ and v^^^. The 
production term is averaged over the planes. Other con- 
tributions in the production term are small due to a 
small corresponding large-scale derivative of the mean 
velocity fleld or they vanish after integration over z [e.g., 

Jq ' C^zUy) dz Ri ii we take into account that the large- 
scale circulation is nearly anti-symmetric relative to its 
central line z — Lz/2 (see Figs. [211S1), where Lz is the 
vertical size of the large-scale circulation] . 

Our measurements of the turbulent heat flux inside 
the coherent structures show that in the regions outside 
the upward and downward large-scale motions the verti- 
cal component of the turbulent heat fiux is very small. 
In particular, the measured vertical turbulent heat flux 
Fz = {ui9) inside the upward and downward large-scale 
motions for AT = 60 K is = 4 K cm s^^, while out- 
side the upward and downward large-scale motions it is 
smaller by at least one order of magnitude. For these 
conditions temperature fluctuations are of the order of 
6 w 1.3 K, the velocity fluctuations are Uz ~ 5.1 cm s~^ 
and the correlation coefficient is 0.6. Turbulence produc- 
tion by the vertical turbulent heat flux Ph = (3 Fz inside 
the upward and downward large-scale motions is of the 
order of Ph ~ 13.3 cm^ (upper estimate). The area of 
the vertical fiow is of the order of 20 % (upper estimate) 
of the area of heat exchangers, which implies that the 
contribution of the turbulent heat flux to the production 
turbulence is 0.2P/i « 2.7 cm^ s~^. 

Measured turbulent energy production due to flow 
shear, averaged over the central part of the flow, is 
P ~ 60.5 cm^ s^^ for AT — 69 K. Since in the inves- 
tigated temperature range the turbulent energy produc- 
tion due to fiow shear P is approximately proportional 
to AT, the magnitude of P at AT = 60 K is of the order 
of P « 50.2 cm^ s~^. Consequently, the contribution of 
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FIG. 5: Dependence of the shear-induced production term 
P = —{ui Uj) V jUi of the turbulent kinetic energy versus the 
dissipation rate e obtained in the experiments with turbulent 
convection in the chamber with A « 2 (triangles for two- 
cell flow pattern and circles for one-cell flow pattern) and 
A « 4 (squares for two-cell flow pattern). The dashed line 
corresponds to P = l.le. The production rate P and the 
dissipation rate e are measured in cm'^ s~^. 

the turbulent heat flux to turbulent energy production 
can be estimated as 5 % (upper estimate). This implies 
that the production of the turbulent kinetic energy in the 
turbulent convection in the regions outside the upward 
and downward large-scale motions is caused by shearing 
motions in the coherent structures. In the vicinity of 
the upward and downward large-scale motions the con- 
tribution to production of the turbulent kinetic energy 
is due to the shear- induced term, P = —(uiUj) V jUi, 
and to the lesser extent by the buoyancy-induced term, 
(3Fz- Our estimations also show that turbulence inside 
the large-scale circulations in thermal convection is pro- 
duced at the expense of a small fraction of thermal energy 
transported by the flow. In particular, these values are 
negligibly small (by a factor of 10~^) in comparison with 
a power transported by thermal convection flow from the 
heater to the cooler. 

In Fig. [5] we show the dependence of the shear- induced 
production term P = —{uiUj) VjUi of the turbulent ki- 
netic energy versus the dissipation rate e in turbulent 
convection. In our experiments we found that P — l.le, 
and the data are scattered by 13 %. This result can be 
explained using a steady state solution of Eq. ^ when 
the production f3 F^ caused by the vertical component of 
the turbulent heat flux is very small. The steady state 
solution of Eq. ^ is P = e. A small deviation from this 
balance equation observed in the experiments can be ex- 
plained by that the formula for the dissipation rate of the 
turbulent kinetic energy includes the empirical constant 
Ce, that can deviate from = 1. Another reason for 
this deviation may be related to the small contribution of 
the turbulent heat flux to turbulent energy production. 

Figure [5] shows the dependence of the shear St^ versus 
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FIG. 6: Measured large-scale shear St^ versus the temper- 
ature difference AT between the bottom and the top walls 
of the chamber obtained in the experiments with turbulent 
convection in the chamber with ^4 ~ 2 (triangles for two-cell 
flow pattern and circles for one-cell flow pattern) and A ~ 4 
(squares for two-cell flow pattern). The temperature differ- 
ence AT is measured in K. 

the temperature difference AT between the bottom and 
the top walls of the chamber, where — iz/u. is the 
turbulent time scale along the vertical direction and S = 
\W zUy\. Here the parameter S is also averaged over the 
probed region in the plane yz. The obtained dependence 
is Stz ~ const (in particular, St. ~ 1 for A w 4 and 
Stz « 0.9 for A ^ 2). This feature is typical for a shear- 
produced turbulence (see, e.g., [6^ [63j). 

In Fig. [7] we show the dependence of the components 
of the measured turbulent velocity u^, Uy and Uz versus 
the temperature difference AT between the bottom and 
the top walls of the chamber. The observed scaling for 
the vertical component of the turbulent velocity is 

Uz (X (AT)0-45, (5) 

while the observed scaling for the horizontal components 
Ux and Uy of the turbulent velocity versus AT has a 
slightly different exponent (w 0.41). In Fig. [5] we show 
the turbulent length scales £x, £y and iz versus the tem- 
perature difference AT. These figures demonstrate that 
the turbulent length scales are nearly independent of the 
temperature difference AT between the bottom and the 
top walls of the chamber. 

It should be also noted that inhomogeneity of the ob- 
served turbulent convection is weak, see Fig. [5] whereby 
the dependencies of the vertical component Uz of the 
measured turbulent velocity and the vertical turbulent 
time Tz = iz/uz versus the horizontal y coordinate are 
shown. 

The scaling ^ can be explained by the following ar- 
guments. Consider the x component of curl of the 
momentum equation for the mean velocity U. In a 
steady-state this equation yields the following estimate 
(g/T,) (VyT) 1/^ A{Vx'U)x. Therefore, the velocity 



7 



OOP 



20 



40 



60 80 AT 



m 



20 



20 



40 



40 



60 80 AT 



60 



c 

80 A7 





□ 


□ 


A 

A 

a 


20 30 40 


50 


60 70 


80 A 


CO 


D 


□ 


A 


20 30 40 


50 


60 70 


80 A 




□ 


□ 


A 

A . 
c 



20 30 40 50 60 70 80 AT 



FIG. 7: Component Ux (panel a), Uy (panel b) and Uz (panel 
c) of the measured turbulent velocity versus the temperature 
difference AT between the bottom and the top walls of the 
chamber (in log-log scale) obtained in the experiments with 
turbulent convection in the chamber with A ~ 2 (triangles 
for two-cell flow pattern and circles for one-cell flow pattern) 
and A « 4 (squares for two-cell flow pattern) . The list-square 
fit for the experimental results is shown by dashed line. The 
turbulent velocity is measured in cm s~^ and the temperature 
difference AT is measured in K. 

shear can be estimated as S" cx {Lz/£)'^ [gT] {VyT)/T^, 
where S ~ VzUy, T* is the characteristic mean temper- 
ature inside the coherent structure and Lz is the ver- 
tical size of the coherent structure. Here we take into 
account that |VzJ7y| :s> \VyUz\, the turbulent viscosity 
ly^ r^u£ = p-jT and ~ SjLl. On the other 

hand, our experiments demonstrate that St = const (see, 
Fig. [6]). Therefore, the turbulent velocity can be esti- 
mated as M (X (3 |Vyr|/r*)i/2^ 

If we assume that |VyT| - {5T)y/Ly and [5T)y cx AT, 
we obtain the following scalings for the turbulent velocity 
field: 

wcx(Ar)i/2, (6) 

and for the dissipation rate of the turbulent kinetic en- 
ergy: 

e = uV2^cx (Ar)3/2. (7) 

In the latter estimate for e we assumed that the turbulent 
length scales £ are nearly independent of the temperature 
difference AT between the bottom and the top walls of 
the chamber (see Fig. [5]). The assumption about the 
linear dependence {5T)y cx AT of the horizontal temper- 
ature difference {5T)y on the vertical temperature differ- 
ence AT, is confirmed by our measurements (see Fig. [TU] 
where the dependence of the horizontal temperature dif- 
ference {ST)y = |VyT|Ly inside one cell of the two-cell 
coherent structure versus the vertical temperature dif- 
ference AT between the bottom and the top walls of 
the chamber with ^ « 4 obtained in our experiments 



FIG. 8: Turbulent length scales (panel a), iy (panel b) and 
£z (panel c) versus the temperature difference AT between 
the bottom and the top walls of the chamber obtained in 
the experiments with turbulent convection in the chamber 
with yl ~ 2 (triangles for two-cell flow pattern and circles 
for one-cell flow pattern) and A ~ 4 (squares for two-cell flow 
pattern). The turbulent length scales are measured in cm and 
the temperature difference AT is measured in K. 




FIG. 9: Dependencies of the vertical component Uz of the 
measured turbulent velocity (upper panel) and the vertical 
turbulent time Tz = £z/uz (lower panel) versus the horizon- 
tal y coordinate obtained in the experiments with turbulent 
convection in the chamber with A ~ 4 (squares correspond to 
measurements at the central yz plane at a; = 15 cm, triangles 
are for the yz plane at a; = 5 cm and circles are for the yz 
plane at a; = 25 cm). The turbulent velocity is measured in 
cm s~^, the turbulent time is measured in seconds and the 
coordinate y in cm. 

is shown). Notably, our experiments also show that the 
shear-induced production of the turbulent kinetic energy 
scales as 

P«ecx(AT)", (8) 
where a = 1.38 for ^ « 4 and a — 1.21 for A w 2 (see 



8 



10 



' □ 



6 
4 
2 



JO 



10 20 30 40 50 60 70 AT" 



P 



90 
70 

50 
40 

30 
20 

10- 



oo. 



20 30 40 50 60 70 80 AT 



FIG. 10: Dependence of the horizontal temperature differ- 
ence {5T)y = \\7yT\Ly inside one cell of the two-cell coherent 
structure versus the vertical temperature difference AT be- 
tween the bottom and the top walls of the chamber obtained 
in the experiments with turbulent convection in the chamber 
with A ~ 4 (the filled squares denote left cell and the un- 
filled squares denote right cell) . Here Ly is the horizontal size 
in the y direction of the one cell of the coherent structure. 
The list-square fit for the experimental results is shown by 
dashed line. The temperature differences {ST)y and AT are 
measured in K. 



Fig. [TT|) . A small difTerence in the exponents in the above 
estimates and the experimental results might be caused 
by the weak dependence of the turbulent length scales 
on the temperature difference AT, and a small deviation 
from the linear dependence {6T)y oc AT. 

In this study we presented the dependencies of the 
turbulent flow parameters versus the temperature dif- 
ference AT between the bottom and the top walls of 
the chamber rather than on global Rayleigh number, Ra 
=f3 (AT) H^/{v k), because of the following reasons. Our 
experiments among the others demonstrate that the life- 
time of large-scale circulations (LSC) is much larger than 
the characteristic turbulence time scale r, and character- 
istic spatial scales of LSC are much larger than the char- 
acteristic turbulence spatial scales. The time of transport 
of a fluid particle through LSC is much larger than the 
time r. This implies that LSC can be described in terms 
of a mean-field approach, i.e., LSC is a mean-field object. 
Therefore, turbulent viscosity and turbulent temperature 
diffusivity, which are much larger than molecular trans- 
port coefficients (the kinematic viscosity, and the tem- 
perature diffusivity, k), become more important. Con- 
sequently, the effective Rayleigh number based on tur- 
bulent viscosity and turbulent temperature diffusivity, is 
an important parameter characterizing the properties of 
LSC. It is plausible to suggest also that when the global 
Rayleigh number is much larger than the threshold value 
required for the excitation of convection, the properties 
of LSC are weakly dependent on the magnitude of the 
global Rayleigh number. 



FIG. 11: Production rate P versus the temperature difference 
AT between the bottom and the top walls of the chamber 
(in log-log scale) obtained in the experiments with turbulent 
convection in the chamber with A ^ 2 (triangles for two-cell 
flow pattern and circles for one-cell flow pattern) and A ~ 4 
(squares for two-cell flow pattern). The list-square fit for the 
experimental results is shown by dashed line. The production 
rate P is measured in cm^ s~^ and the temperature difference 
AT is measured in K. 

We show the dependencies of the turbulent flow pa- 
rameters versus the temperature difference AT between 
the bottom and the top walls of the chamber with two 
aspect ratios A 2 and A ^ 4 in the same figure. When 
the temperature difference AT is varied, the kinematic 
viscosity v{T)^ the temperature diffusivity k(T) and the 
parameter /?(T) which depend on the mean fluid tem- 
perature, are also changed. Therefore, several variables 
(AT, u{T), k{T) and /?(T)) in the expression for the 
global Rayleigh number are changed when the temper- 
ature difference AT is varied. Another reason for using 
the temperature difference AT rather than the global 
Rayleigh number as independent variable is that we de- 
termined theoretically the scalings of the production rate 
P and turbulent (rms) velocity u with the temperature 
difference AT rather than on global Rayleigh number, 
Ra. However, for comparing our results with other stud- 
ies we also presented some dependencies of the turbu- 
lent flow parameters versus global Rayleigh number (see 
Fig. [HI). 



IV. DISCUSSION AND CONCLUSIONS 

We investigated the effect of large-scale coherent struc- 
tures on global properties of turbulent convection in lab- 
oratory experiments. The turbulent convection for the 
Rayleigh numbers varying from 5 x 10^ to 10®, was stud- 
ied in air flow in a rectangular chamber with aspect ra- 
tios A K, 2 and A « 4. The large-scale coherent struc- 
tures comprise the one-cell and two-cell flow patterns. It 
was demonstrated that turbulent convection with large- 
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FIG. 12: Dependencies of the vertical component Uz of the 
measured turbulent velocity (upper panel) and the production 
rate P (lower panel) versus the global Rayleigh number Ra 
(in log-log scale) obtained in the experiments with turbulent 
convection in the chamber with A ^ 2 (triangles for two-cell 
flow pattern and circles for one-cell flow pattern) and A ~ 4 
(squares for two-cell flow pattern). The list-square flt for the 
experimental results is shown by dashed line for A « 2 and 
by dashed-dotted line for A « 4. The production rate P is 
measured in cm^ s~'^ and the turbulent time is measured in 



ment with theoretical predictions. 

Equation ([1]) that is valid for a a buoyancy-produced 
turbulence, yields much smaller values for the turbulent 
kinetic energy than those observed in our experiments. 
The production of the turbulent kinetic energy in tur- 
bulent convection, caused by the buoyancy flux, even in 
the upward and downward flow regions is much smaller 
than the production rate due to the shear motions inside 
the large-scale coherent structures. Inside the coherent 
structures where the vertical heat flux is very small, the 
production of the turbulent kinetic energy, caused by the 
buoyancy flux, is by several orders of magnitude smaller 
than the production rate due to the shear motions. 

The large-scale coherent structures are formed due to 
large-scale instability caused by a redistribution of the 
turbulent heat flux due to non-uniform motions. As soon 
as the large-scale coherent structures are formed, there 
is a strong influence of the shear produced by the coher- 
ent structures. In particular, the main contribution to 
production of turbulence is due to shearing motions in- 
side the large-scale coherent structures. The magnitude 
of the large-scale shear is determined by the horizontal 
mean temperature gradient inside the large-scale circu- 
lation. These features are different from those which are 
typical for the classical laminar convection. 



scale coherent structures is a shear-produced turbulence. 
We determined the dependence of global parameters (the 
production and dissipation of turbulent kinetic energy, 
the turbulent velocity and integral turbulent scale, the 
large-scale shear, etc.) of turbulent convection on the 
temperature difference between the bottom and the top 
walls of the chamber. These dependencies are in an agree- 
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